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Abstract 



We consider the problem of approximating the partition function of the hard-core 
model on planar graphs of degree at most 4. We show that when the activity A is suffi- 
ciently large, there is no fully polynomial randomised approximation scheme for evaluating 
the partition function unless NP = RP. The result extends to a nearby region of the pa- 
rameter space in a more general two-state spin system with three parameters. We also give 
a polynomial-time randomised approximation scheme for the logarithm of the partition 
function. 

1 Introduction 

A spin system is a model of particle interaction on a graph. Every vertex of the graph is 
assigned a state, called a spin. A configuration assigns a spin to every vertex, and the weight 
of the configuration is determined by interactions of neighbouring spins. 

In this paper, we consider the following two-spin model, which applies to spin systems 
on a graph G = (V,E). The model has three parameters, f3, 7 and A. It is easiest to 
view these as non-negative rationals for now — we will be slightly more general later. A 
configuration a: V{G) — > {0, 1} is an assignment of the two spins "0" and "1" to the vertices 
in V. The configuration a has a weight wq{o\ which depends upon /3, 7 and A. Let b(o~) 
denote the number of edges (u, v) of G with a(u) = cr(v) = 0, let c(<r) be the number of 
edges (u, v) of G with a(u) = a{y) = 1 and let t{a) be the number of vertices u of G with 
o~(u) = 1. Then wg(o~) = ft b ^ 7 C ( CT ) A^ CT ) . The partition function of the model is given by 



Yl w g(o-). 

<t:V(GH{0,1} 



Two important special cases are 



• the case f3 = 1, 7 = 0, which is the hard-core model, and 



• the case /3 = 7, which is the Ising model 
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The hard-core model [3] is a model of a gas in which vertices are either occupied by a 
particle (in which case they have spin 1) or unoccupied (in which case they have spin 0). 
The particles cannot overlap and adjacent vertices are close together, hence 7 = 0. The Ising 
model is a model of ferromagnetism. In this paper we study the hard-core model and a region 
of nearby two-state spin systems. 

1.1 Previous work 

Evaluating Zp^ t \{G) is a trivial computational problem if [3j = 1, because the partition 
function factors. In other cases, the complexity of evaluation has been studied in detail. 
When A = 1, the problem of computing the partition function on planar A-regular graphs 
is called P1-Ho1a(o,&) in [7], where a corresponds to f3 and b corresponds to 7. There is a 
dichotomy [3 Theorem 1]: for non-negative a,b, the problem P1-Ho1a(o, b) can be computed 
in polynomial time in the trivial cases oh = 1 and a = b = 0, and in the case of the Ising 
model with no external field, a = b. In all other cases the problem is #P-hard. 

A standard transformation extends this dichotomy to arbitrary A > 0. Consider a con- 
figuration a: V{G) — > {0, 1} of a planar A-regular graph G. Counting the number of edges 
adjacent to a "1" spin in two ways, we have A£(a) = 2c(a) + (\E(G)\ — b(a) — c(a)). Therefore, 

which is as hard to compute as P1-Ho1a(/3A _1//a , jX 1 ^). Suppose /3 and 7 are not both 0. 
Unless A = 1, we have either /3A" 1 / 3 / 7A 1 / 3 or /3A" 1 / 4 / 7A 1 / 4 . If £7 ^ 1 then in either 
case, we can conclude from above that evaluating Zp-y \{G) is #P-hard when the input G is 
restricted to be a planar graph of degree at most 4. 

Since the complexity of exactly evaluating the partition function is intractable, much effort 
has focussed on the difficulty of approximately evaluating the partition function for a given 
set of parameters j3, 7 and A. 

The complexity of approximating the partition function of the hard-core model and the 
Ising model in general (not necessarily planar) graphs is well-understood. The Gibbs measure 
is the distribution on configurations a: V(G) {0,1} in which the probability of configu- 
ration a is proportional to wg(o~). This notion of Gibbs measure extends to certain infinite 
graphs, for example infinite regular trees, where it may or may not be unique. For the hard- 
core model, there is a critical point A C (A) = (A - 1) A -V(A - 2) A such that the infinite 
A-regular tree has a unique Gibbs measure if and only if A < A c . An important result of 
Weitz [26j showed that, in every graph with maximum degree at most A, the correlations 
between spins in the hard-core model decay rapidly with distance as long as A < A c . As a 
result, he gives [26, Corollary 2.8] a fully-polynomial (deterministic) approximation scheme 
(FPTAS) for evaluating the hard-core partition function on graphs of degree at most A for 
any A < A c . By contrast, Sly and Sun [241 Theorem 1] (see also the earlier hardness results of 
Sly [25] and Galanis et al. [10]) show that, unless NP = RP, there is no fully-polynomial ran- 
domised approximation scheme (FPRAS) on A-regular graphs (for A > 3) for any A > A C (A). 
Thus, the difficulty of approximation is resolved, apart from at the boundary A = A C (A). 

We say that the two-spin model is ferromagnetic if ^7 > 1 and antiferromagnetic if 
/?7 < 1. For the antiferromagnetic Ising model, Sinclair et al. [231 Corollary 1.2] show that 
there is a FPTAS for evaluating the Ising partition function on graphs of degree at most A 
for any choice of parameters (3 and A which is in the interior of the uniqueness region of the 
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A-ary tree. By contrast, Sly and Sun |24|, Theorem 2] show that, unless NP = RP, there is no 
FPRAS on A-regular graphs (for A > 3) if /3 and A are outside the uniqueness region. (So, 
once again, the situation is fully resolved, apart from the boundary.) The result of Sinclair et 
al. |23|, Corollary 1.3] extends to general anti- ferromagnetic two-state spin systems in regular 
graphs, and also in a somewhat wider class of graphs. 

For general anti-ferromagnetic two-state spin systems, the best positive result that is 
known is due to Li, Lu, and Yin [19]. They use a stronger notion of correlation decay than 
Weitz, which enables them to obtain a PTAS, even for graphs with unbounded degree. They 
show [T9l Theorem 2] that for any finite A > 3, or for A = oo, there is an FPTAS for the 
partition function of the two-state spin system on graphs of maximum degree at most A if the 
parameters of the system are antiferromagnetic, and for every d < A, they lie in the interior 
of the uniqueness region of the infinite d-regular tree. By contrast [T9l Theorem 3], the results 
of Sly and Sun imply that, for any finite A > 3, or for A = oo, unless NP = RP, there is 
no FPRAS for the partition function of the two-state spin system on graphs of maximum 
degree at most A if the parameters of the system are antiferromagnetic, and for some d < A, 
they lie outside the interior of the uniqueness region of the infinite cf-regular tree. Thus, 
the approximation complexity is resolved in the antiferromagnetic case, apart from at the 
boundaries of the uniqueness regions. Note that the result of Sun and Sly was independently 
discovered by Galanis, Stefankovic and Vigoda [11] for the case A = 1. 

The situation is not completely resolved in the ferromagnetic case. Building on Jerrum and 
Sinclair's FPRAS for the ferromagnetic Ising model [16], Goldberg, Jerrum and Paterson [15] 
gave an FPRAS for the ferromagnetic two-spin model which applies if /3 > 7 and A < 
\[Wp~) (or, equivalently, if (3 < 7 and A > \fWfl)- The approximation applies without these 
constraints on the parameters if the input is a regular graph. 

For the hard-core model, an important issue which arises in statistical physics is approx- 
imating the partition function for planar graphs, including regular lattices. While (as far as 
we know) there were no hardness results for this problem (until this paper) the complexity of 
particular algorithms have been studied. For example, Randall [21] showed that a particular 
MCMC algorithm provides a bad approximation on subsets of Z, 2 , because Glauber dynamics 
mixes slowly when A > 8.066. (By contrast, results of Restrepo et al. [22] showed that the 
mixing time is 0(n log n) when A < 2.3883, and that Weitz's algorithm |26j gives a (deter- 
ministic) fully-polynomial-approximation scheme in this case.) Recently tree decompositions 
of planar graphs have been used to give FPTASes for certain partition functions on planar 
graphs — see [27] . 

1.2 Our contribution 

Our objective is to determine whether approximating the partition function of the hard-core 
model is computationally intractable on planar graphs for sufficiently large A. It turns out that 
this is so. Our main result (see Theorem [1]) is that, for a wide range of two-spin parameters, 
there is no FPRAS, even for planar graphs with degree at most 4. The applicable range of 
parameters includes the hard-core model with A > 312. Thus, we show that approximation 
is difficult for this problem (see Corollary [2]) . 

An interesting difference between the general case and the planar case is that, in general, 
it is difficult to approximate the logarithm of the partition function, a quantity which has 
physical significance and is called the mean free energy. Sly and Sun (see the proofs of 
Theorems 1 and 2 in |24J) showed that there is a fixed c > 1 such that no algorithm can 
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approximate Zp~\{G) within a factor c^ v ^ G ^ unless NP=RP. By contrast, we show (see 
Theorem [3J that, in the planar case, there is a polynomial-time approximation scheme for 
logZ ( g 7) ^(G) (which implies an approximation within a factor c^ v ^ G ^ since Z^ 7 ^(G) is at 
most C\ V ( G ^ for a quantity C which depends only on (3, 7 and A). 

At a high level, our hardness result is a reduction from the optimisation problem of 
computing a maximum independent set in a cubic planar graph G to the problem of estimating 
the partition function of a much larger graph, which is constructed from G. Each vertex of G 
is represented by a gadget which is a "wrapped" rectangular lattice C v (see Figure [I]) . Similar 
to previous results of Goldberg and Jerrum [13], and Sly [25], and Sly and Sun [23], we exploit 
the phase transition of the gadget to enable a reduction from a hard optimisation problem. 
The optimisation problem from which we start (computing a maximum independent set in a 
cubic planar graph) plays a similar role to that of the maximum cut problem in the reduction 
of Sly and Sun [24J. However, there is a key difference. Since, as we discuss below, it turns 
out that the logarithm of our partition function is efficiently approximable, it is therefore 
necessary that the optimisation problem from which we start is also easy to approximate 
(otherwise, we would get a contradiction). This means that our reduction has to be more 
carefully tuned — the approximation of the partition function has to allow us to exactly solve 
the optimisation problem. 

A key technical challenge in the proof is to characterise the Gibbs distribution of the 
two-spin model on the lattice gadget. We show that the spins of the vertices do exhibit long- 
range correlation. In fact, the gadget is almost always in one of two phases. Each of these 
phases are equally likely. Also, conditioned on the phase, the spins of certain vertices along 
the boundary of the gadget are nearly independent, and their distribution can be determined. 
Thus, although there is long-range correlation between spins, all of the correlation is captured 
by the phase. Conditioned on the phase, the spins are not very correlated. The analysis of 
the Gibbs distribution of the gadget uses contour arguments adapted from Dobrushin [8j and 
Borgs et al. [5j. Randall's slow- mixing result is also based on contour arguments. 

In statistical physics it is sometimes useful to approximate the logarithm of the partition 
function, even when the partition function itself cannot be approximated (for example, in the 
situation of Theorem [1]). Bandyopadhyay and Gamarnik [2 J have shown how to estimate the 
logarithm of the partition function of the hard-core model when A is small and the graph is 
regular, with large girth. They show that, in this case, the approximate value does not depend 
on the graph, given its degree and size! We give (Theorem [3|) an approximation scheme for 
the logarithm of the partition function which applies to all planar graphs, for sufficiently 
large A. The algorithm is based on the decomposition technique that Baker [1] used to give 
approximation schemes for optimisation problems on planar graphs. There is a parameter k 
which is governed by desired approximation quality. The graph G is decomposed into pieces 
which are /c-outerplanar, and therefore have bounded tree-width. The partition functions of 
these pieces can be calculated directly using an algorithm of Yin and Zhang [27] . These are 
combined to give the estimate. 

2 Preliminaries and statement of results 

In our main result, we will assume that the parameters /?, 7 and A satisfy the following 
conditions. 

A > 1, P > 1 > 7 > 0, 07 < 1, /3A" 1/4 < 0.238, and 7A 3/8 < 0.238. (1) 
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Note that these conditions are satisfied by the hard-core model when A > 312 (by setting 
P = 1 and 7 = 0). 

The notion of a fully polynomial randomised approximation scheme (FPRAS) is defined 
in Section [3j Following [13^ 114] , we say that a real number z is efficiently approximate if 
there is an FPRAS for the problem of computing z. For fixed efficiently approximable reals 
(3, 7 and A satisfying ([1]), we consider the problem of (approximately) computing Z^^ t \(G), 
given an input graph G. In order to make our (negative) result as strong as possible, we 
restrict the input G to have degree at most 4 as well as being planar. Thus, we study the 
following computational problem. 

Name DegreeFourPlanarTwoSpin(/3, 7, A). 

Instance A planar graph G with maximum degree at most 4. 

Output The value Zp n< \{G). 

Our main result is the following. 

Theorem 1. Suppose that f3, 7 and A are efficiently approximable reals satisfying ([T|). There 
is no FPRAS for DegreeFourPlanarTwoSpin(/3, 7, A) unless NP C BPP. 

The inclusion NP C BPP would imply NP = RP [18l Theorem 2]. So, for any fixed /3, 7 
and A satisfying (JTJ), there is no FPRAS for DegreeFourPlanarTwoSpin(/3, 7, A) unless 
NP = RP. 

Of course, our result has an immediate consequence for the problem of approximating 
the partition function in the hard-core model. Thus, Theorem Q] implies Corollary [2] for the 
following computational problem. 

Name DegreeFourPlan arHardCore( A) . 

Instance A planar graph G with maximum degree at most 4. 

Output The value Zi t o t \(G). 

Corollary 2. Suppose that A > 312 is an efficiently approximable real. There is no FPRAS 
for DegreeFourPlanarHardCore(A) unless NP C BPP. 

Despite Theorem [H we show that the logarithm of the partition function can be approxi- 
mated. In particular, we study the following computational problem, where, for concreteness, 
we use the natural logarithm (to the base e). 

Name PlanarLogTwoSpin(/3, 7, A). 
Instance A planar graph G. 
Output The value log(Z / g >7 : \(G)). 

Our result is that there is a polynomial-time randomised approximation scheme (PRAS) 
for PlanarLogTwoSpin(/3, 7, A). A polynomial-time randomised approximation scheme is 
a more liberal notion of approximation than the notion of an FPRAS. See Section [3] for a 
definition. 
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Theorem 3. Suppose that (3, 7 and A are efficiently approximable reals satisfying /3 > 1 > 
7 > and A > 1. There is a PRAS for PlanarLogTwoSpin(/3, 7, A). 

The randomness used by the algorithm promised by Theorem [3] is only needed to ap- 
proximate the parameters /3, 7 and A. If these are deterministically approximable, then the 
approximation is deterministic. 

We will need some notation to refer to the Gibbs distribution of the two-spin model on 
a graph G, which is the distribution in which the the probability of each configuration is 
proportional to its weight. We will use ctq to denote a random configuration drawn from this 
distribution. Thus, for any configuration a: V(G) — > {0,1}, 

Pr(er G = a) = wa{cr)/Zp at x(G). 

(In general, as here, we use boldface for the random variable and normal type for the values 
that it takes on.) Finally, given a subset S of V(G) and a configuration a : V(G) — > {0, 1}, 
let a(S) : S — > {0, 1} denote the configuration induced by a on S. 

3 Polynomial Randomised Approximation Schemes 

Most of this section is taken from [14] and can be skipped by readers who are already familiar 
with randomised approximation schemes. 

A randomised approximation scheme is an algorithm for approximately computing the 
value of a function / : S* — >• R. (Here, £ is a finite alphabet, and inputs to / are represented as 
strings over this alphabet.) The approximation scheme has a parameter e > which specifies 
the error tolerance. A randomised approximation scheme for / is a randomised algorithm that 
takes as input an instance x G S* (e.g., for the problem DegreeFourPlanarTwoSpin(/3, 
7, A), the input would be an encoding of a planar graph G) and a rational error tolerance 
£ € (0, 1), and outputs a rational number z (a random variable of the "coin tosses" made by 
the algorithm) such that, for every instance x, 

Pr [e- £ f(x) < z < e e f(x)] > j . (2) 

The randomised approximation scheme is said to be a polynomial randomised approximation 
scheme or PRAS if, for each e, its running time is bounded by a polynomial in \x\. It is said 
to be a fully polynomial randomised approximation scheme, or FPRAS, if its running time is 
bounded by a polynomial in |x| and e . 

Note that the quantity | in Equation (|2|) could be changed to any value in the open interval 
(^, 1) without changing the set of problems that have randomised approximation schemes [TTl 
Lemma 6.1]. In fact, in the proof of Theorem[Tl we will assume that our FPRASes have failure 
probability at most 1/15. 

The notion of an FPRAS is a particularly robust notion of approximability for parti- 
tion functions. For such approximations, the existence of a polynomial-time algorithm that 
achieves a constant-factor approximation actually implies the existence of an FPRAS. The 
same argument that we gave to illustrate this point for the Potts model [14] also applies 
to the setting of this paper. For any graph G, denote by k ■ G the graph composed of k 
disjoint copies of G. Then Zp a ^\(k ■ G) = Zp a ^\{G) k . So, setting k = 0(e _1 ), a con- 
stant factor approximation to Zp~ u \(k ■ G) will yield (by taking the A;th root) an FPRAS for 
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DegreeFourPlanarTwoSpin(/3, 7, A). Clearly, an approximation within a polynomial fac- 
tor would also suffice. Note that the same argument does not necessarily apply to log-partition 
functions. 

4 The Gadget 

We will assume throughout this section that /3, 7 and A satisfy ^TJ, so we do not keep repeating 
this condition in the statement of our lemmas. 

The gadget C u has vertex set V(C U ) = Z/2vZ x {0, . . . , u}. Vertices (x, y) and (x 1 , y') are 
adjacent in C v if 

• y = y' and x = x' ± 1 (where of course, the arithmetic is modulo 2v since x and x' are 
in Z/2i/Z), or 

• x = x' and y = y' ± 1. 

Let E(C U ) denote the set of edges of C v . See the leftmost picture in Figure [TJ 

(5,0) (5,0) 




(0,0) (0,0) 

Figure 1: C5, and the vertex subsets #1,2, and i?o,5- 
4.1 Goalposts and keyholes 

Given a vertex (x, 0) G V(C U ) and a value m G {0, . . . , v} let -B Xjm be the set containing the 
vertices on the rectangular (goalpost-shaped) path at distance m around the terminal. In 
particular, let 

B x ,m= (J {(x -m,j),(x - j,m),(x +j,m),(x + m,j)}. 

0<j<m 

Again, the arithmetic is done modulo 2v since x G Z/2^Z. See the middle picture in Figure[TJ 
When m = v, the vertices in {(x — m,j) | < j < m} coincide with the vertices in 

{{x + m,j) I < j < m} so -B^m becomes the "keyhole" which is depicted in the rightmost 

picture of Figure [U (for x = 0). 

We shall often be working with configurations on gadgets. For convenience the notation 

<jc v will be contracted to cr u , and no confusion should result. 
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4.2 Parity-0 ones and parity-1 ones 

We say that a vertex (x, y) G V{C V ) has parity 1 if x + y is odd, and that it has parity 
otherwise. Suppose that S is a subset of V(C V ) and that s G {0, 1}. We say that a(S) has 
parity-s ones if {(x,y) G S | a(x,y) = 1} is exactly the set of parity-s vertices in S. 

4.3 Idealised probabilities 

Define 

p~ = limsupPr(<7v(0, 0) = 1 | ct u (Bq^) has parity-0 ones), and 

v— >oo 

p^ = limsupPr(<T I/ (l,0) = 1 | a v {B\ tV ) has parity-0 ones). 

v— >oo 

The notation p = is meant to connote that we are looking at the probability of a 1 at a 
vertex of parity s, conditioned on certain parity-t ones, where s = t; for p^ we are interested 
in s ^ t. As we shall see later, it will turn out that p = > p^. This is a non-trivial fact 
about the spin system: if there were no long-range correlations, we would have p = = . The 
following straightforward lemma is also useful. 

Lemma 4. p^ > and p = < 1. 

Proof. Suppose v> 2. Consider vertex (1, 0) of C v . Let S = {(1,0), (2,0), (0,0), (1, 1)} be the 
set containing (1, 0) and its immediate neighbours. Let S' = {(—1, 0), (0, 1), (1, 2), (2, 1), (3, 0)} 
be the set containing the neighbours of S. Given any a : S' — > {0, 1}, 

Pr(ov(,S) has parity-1 ones | a u (S r ) = a) > A/16A 4 /3 10 > 0. 

Now let S" = {(-1, 0), (0, 1), (1, 0)} be the neighbours of (0, 0). Given any a : S" -> {0, 1}, 

PrK(0,0) = 1 | a u (S") = a)< -A- < 1. 

1 + A 

□ 

The events that a u (Bo jU ) has parity-0 ones and that <r v {B\ iV ) has parity-0 ones have low 
probability, so it may seem strange to condition on these events, but the purpose of this 
conditioning is to identify two phases of the idealised gadget. We will refer to certain vertices 
(x, 0) of C u as "terminals", and it will turn out to be the case that the spins of these terminals 
are nearly independent of each other in the distribution of cr u . 

We will study the distribution that a u induces on the terminals by considering an idealised 
distribution with two phases. In each of these two gadget phases, the spins of the terminals 
will be chosen independently. Some terminals will be assigned spin 1 with probability p = and 
others will be assigned spin 1 with probability p^ . This will be explained further in the next 
section. 

4.4 Terminals 

Fix positive integers d and k. Let v = 2dk and let denote the gadget C u . We will work 
with Ck,d for the rest of the paper. We will use both notations, C v and Ck,d, depending on 
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whether we want to emphasize the role of v or the role of k and d. Similarly, the alternative 
notations, (Tk,d and cr u will be used as convenient. 

Some of the vertices around the boundary of Ck,d (2k of them) are designated as "termi- 
nals". The set of "parity- 1 terminals" is 

Tl d = {(4jd + l,0)\0<j<k-l}. 

The set of "parity-0 terminals" is 

Tl d = {(4jd + 2d,0) \ 0<j<k-l}. 

Let Tk^d = T k d UT k d denote the set of terminals. 

For parity s £ {0, 1}, let [i s k d be the distribution on configurations a : Tk t d — > {0, 1} in 
which the spin of each terminal is chosen independently as follows: For each parity-s terminal 
(x,0), set a(x,0) = 1 with probability p = (and set a(x,0) = otherwise). For each terminal 
(x,0) with parity 1 © s, set a(x, 0) = 1 with probability p^ (and set a(x,0) = otherwise). 

Informally, the distribution /x? , will be relevant when an idealised gadget is in a phase 
which prefers 1-spins at parity-s terminals. In this distribution, the probability that a terminal 
is given spin 1 is higher if the terminal has parity s than if it has parity 1 © s. 

Let fik,d be the distribution on configurations a : T^ : d —> {0, 1} given by /J>k,d(o-) = 
(/ i fcd( <T ) + Mfe d( <7 ))/2- We will show that, provided that d is sufficiently large, the distri- 
bution of (Tk,d{Tk,d) is close to 

Thus, the gadget can be thought of informally as having two phases, phases and 1. We 
will show that the gadget almost always occupies one of these two phases, and they occur 
with equal probability. In phase 0, the distribution of &k,d^k,d) 1S close to d . In phase 1, 
the distribution of o-k,d{Tk,d) is close to n\ d . 

Proposition 5. There is a c > 1 such that, if d is a sufficiently large multiple of 16, k is an 
integer greater than or equal to 1 and r is a configuration r : Tk t d —> {0, 1}, then 

| Pr(<r M (T M ) = r) - /i m (t)| < c~ d k 2 . 

Proposition [5] is established at the end of this section. We will use contour arguments 
adapted from Dobrushin [S] and Borgs et al. [5j. The outline of the argument is as follows. 
We first define "contours" in Section 14.51 We show, in Section 14.61 that long contours are 
unlikely. In Section T4.71 we show that, in the absence of long contours, the spins of terminals 
are nearly independent. With high probability, the gadget has a phase s and there is a 
boundary around each terminal, whose spins are consistent with s. Conditioned on s, the 
distribution of the spins of the terminals is close to fi s k d . 

4.5 The Dual Gadget, trails, and contours 

The dual gadget C* has vertex set V(C*) = {x + | | x G Z/2vZ} x{y + \\y e {-1, . . . , v}}. 
Vertices (x,y) and (x',y') are adjacent in C* if 

• y = y' and y ^ {— i, ^+^}, and x = x'±l (where of course, the arithmetic is modulo 2u), 
or 

• x = x' and y = y' ± 1. 
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Figure 2: Part of C3 and C3; solid lines are edges of C3, dashed lines are edges of C3. The 
red thickened lines are a dual pair of edges. 

E(C*) denotes the edge set of C*. This is illustrated in Figure [3 

There is a bijection called "duality" between edges of C v and edges of C*. In particular, 
the dual of edge e = ((x,y),(x + l,y)) of C v is e* = ((x + -%,y — |), (x + \,y + 5)) and 
the dual of edge e* is e. Similarly, the dual of edge / = (x, y + 1)) of C u is /* = 

((x — ^,y + ^),(x + ^,y + ^)) and the dual of /* is /. We use the * operation to move between 
an edge and its dual, so every edge e satisfies (e*)* = e. 

A trail in C* is a sequence g = v±, . . . , Vj of vertices in V(C*) such that each pair (vi, i>i+i) 
is an edge of C*, and no edge is used twice. A contour is a trail g = V\, . . . ,Vj in C* satisfying 
one of the following: 

• vi = Vj, or 

• The y-coordinate of v\ and the y-coordinate of Vj are both in {— |, v + |}- 

The length of g is j — 1. We say that 5 is a cross contour if the y-coordinate of v\ is — ^ and 
the y-coordinate of Vj is v + ^ (or vice- versa). A cross contour goes from one boundary of 
the gadget to the other. We say that every other contour is a simple contour. 

Given a : V(C U ) — > {0, 1}, let a* be the set of edges of C* which are dual to monochromatic 
edges. In particular, a* = {(u,v)* G E(C*) \ a{u) = o(v)}. 

Definition 6. A contour of a is a contour g = v\, . . . ,Vj satisfying the following two proper- 
ties. 
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• The edges of g are monochromatic: That is, for all 1 < i < j, (i^Uj+i) € cr* . 

• The contour g always turns at degree-A vertices: That is, for all 1 < i < j, if four edges 
of a* meet at vertex Vi, then Vi—i and Vi+i differ in both the x component and the y 
component. Similarly, if four edges of a* meet at v\ = Vj then v<i and Vj-i differ in 
both the x component and the y component. 

Note that contours of a cannot cross each other, though two contours can share a vertex 
without crossing. Also, two contours can have a common portion (before turning off in two 
different directions). Finally, every edge of a* is contained in at least one contour of a. 

Let a : V(C U ) — > {0, 1} be a configuration, and let g be a contour of a. We say that a 
vertex u € V(C U ) is adjacent to g if there is an edge (u, v) £ E(C V ) such that e* € g. The set 
of vertices adjacent to g can be written as the union of two sets, L{g) and R(g), where L(g) 
is the set of vertices which are on the left (relative to the direction of travel) when we follow 
the trail g from v\ to Vj, and R{g) is the set of vertices which are on the right (relative to the 
direction of travel) . See Figure [3l 























L 












L - 


T-7*77T 

liaiii 


- L 








L - 


:::: > ::: 

::: r :: 




H 


-> 


- L 








L - 


jiijii 
::: < :::: 


llllRl 

y.'.'.y.: 




- L 










L 


L 





















Figure 3: Left and right vertices of a contour of a. The shaded squares represent vertices of 
C„ with parity- 1 ones and the unshaded squares represent vertices of C u with parity-0 ones. 

A key property of contours is the following. 

Lemma 7. Let a : V(C U ) — > {0, 1} be a configuration, and let g be a contour of a. Then for 
some s 6 {0, 1}, a(L(g)) has parity-s ones and a(R(g)) has parity-(l(&s) ones. 

Proof. Pick s € {0, 1} such that the vertex on the left as we go from v% to x>2 has parity-s ones. 
By induction on i, we will show that for each i the vertex on the left as we go from Vi to Vi+i 
has parity-s ones. Suppose without loss of generality that the edge increases the x- 

component (the other three cases are similar). So i>j_i = (x — h,y + ^) and V{ = (x + + 
Since g is a contour of a, a(x, y) = a(x, 7/ + l) = s©x0y. There are three cases. 

1. Wj+i = (x + |, y + |). In this case the vertex (x, y + 1) is still on the left as we go from 

Vi to V i+ i. 

2. Vi + i = (x + — i). In this case the vertex (x + l,y) is on the left as we go from V{ 
to Vi+i, but since g is a contour of a we have a(x + l,y) = o~(x,y) = cr(x,y + 1). So 
(x + l,y) has parity-s ones. 



11 




X X + 1 X X + 1 X X + 1 



Figure 4: Cases 1, 2, and 3. Black squares have the same spin as a(x, y); white squares have 
the opposite spin, and hatched squares can be either. 

3. Vi+i = y+|). In this case (x + l,y + l) is on the left, and (x + 1, y) is on the right. 

Since g is a contour of a, a(x + l, y) = a(x + 1, y + 1), and we know a(x, y) = a(x,y + l). 
Since the contour did not turn, the vertex (x + \ ,y + \) cannot have degree 4 in a*, so 
a(x + 1, y + 1) = s © x © y © 1, so (x + 1, y + 1) has parity-s ones. □ 

The following lemma allows wa{o~) to be expressed more easily in terms of the contours 
of a. Suppose v > 2. A side vertex of C v is a vertex (x, y) € V(Cj,) with y = or y = v. A 
side edge is an edge in E{C V ) between two side vertices. 

Lemma 8. Fix u > 2 and a configuration a : V{C V ) — > {0, 1}. Xei &'(<t) be the number of 
side edges (u,v) of C v with a(u) = a(v) = and let c'(cr) be the number of side edges (u,v) 
of C v with a(u) = o~{v) = 1. Then 

£(a) = \(c(a) - b(a)) + |(c'(a) - b'(a)) + u(y + 1). 

Proof. Let ^'(tr) be the number of side vertices u with o~(u) = 1. Let £" be the set of all side 
edges in C^. By double-counting pairs (u, (u,v)) with a(u) = 1 and (u,v) £ E{C V ), 

(\E(C„)\ - b(a) - c(a)) + 2c(a) = Ai{a) - £'(a). 

By double-counting pairs (u, (u,v)) with cr(u) = 1 and (u,v) € E', we have 

- b'(a) - c'{a)) + 2c' {a) = 2l\o). 

Rearranging gives £{a) = \{c(a) - b(a)) + §(c'(cr) - + \\ E {Cv)\ + \\E'\. Consider the 

configuration with alternating 0s and Is given by a(x,y) = x © y. For this configuration a, 
we have b(a) = c(a) = b'(a) = c'(a) = and £(a) = v(y + 1), so the constant term 
\\E(C V )\ + \\E'\ is 1/(1/ + 1). ' n 

4.6 Long contours are unlikely 

Lemma 9. There is a c > 1 such that, for all sufficiently large h, all u > 2, and all U C 
V(C* U ), 

Pr(cr I/ has a simple contour of length at least h starting in U) < \U\ c~ h . 

Proof. Suppose that g is a simple length-r contour of a configuration a: V(C V ) — > {0,1}. 
Consider the connected components of the graph [y(C v ),E(C v ) \ {e* | e G 5}). We say 
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that a component is "left" if it contains at least one vertex in L(g) (but no vertices in R(g)). 
We say that it is "right" if it contains at least one vertex in R(g) (but no vertices in L{g)). 
Every component is either left or right. Let S be the set of vertices in left components. Let 
S = V(C V ) \ S. Let S' = {(x, y) G S \ (x — 1, y) G S}, where, as usual, the arithmetic on x is 
done modulo 2u. 

Suppose that a(R(g)) has parity-s ones. By Lemma [71 this is true for some s G {0,1}. 
Define a configuration a 9 : V{C U ) — > {0,1} as follows: a 9 (S) = a(S), cr 9 (S') has parity-s 
ones, and, for every (x, y) G S \ S', a 9 (x, y) = a(x — 1, y). 

Note that a i— > a 9 is a map from the set of configurations a with g as a contour to the set 
of all configurations; further, it does not lose information, and hence is injective. Note also 
that (a 9 )* is the same as a*, but with g removed and with the contours in S shifted by one. 
By Lemma El 

WCv {a) = (/3A- 1 / 4 )^( 7 A 1 / 4 ) c ^A( c '^- fe '^)/ 8 U ;c,(^) 

where b(g) , c(g) , b' (g) , d \g) are the contributions to b(a) , c(a) , b' (a) , d (a) coming from edges 
whose duals are in g. As the map a i-> a 9 is injective, 

Pr(<? C a* v ) < ( / 3A- 1 /4)%)( 7A V4 )C (< ? ) A ( C '( 9 )-;>'( 9 ))/8 

<(^A- 1 /4)%)( 7A 3/8 )C ( S ) ; 

where we have used the facts that A > 1 and d{g) < c{g). There are at most \U\ 3 r relevant 
contours of length r in total {\U\ choices of starting point, and at most three different directions 
at each step), so 

Pr(cr i/ has a simple contour of length at least h starting in U) 

< |[/|^3 r max(/3A- 1/4 , 7 A 3/8 ) r 

r>h 

< , (3 x 0.238) fe 
~ ' ' 1 - 3 x 0.238' 



There is a c > 1 such that j^g ^ 3 ^ < c h for all sufficiently large h. □ 

Lemma 10. Let i G {0, 1}. For every v > 2, and every simple contour g of length r, 

Pr(<7 is a contour of cr u \ rT v {B^ v ) has parity-^ ones) < max(/3A -1 ^ 4 , 7A 3 ^ 8 ) r . (3) 

Furthermore, there is a c > 1 such that, for all sufficiently large h, all v > 2, and all 
U C V(C U ), the conditional probability that <j u has a simple contour of length at least h which 
contains an edge whose dual connects two vertices in U, conditioned on the fact that cr v (Bi jV ) 
has parity-0 ones, is at most \U\ c~ h . 

Proof. The proof of ^ is similar to the first half of the proof of Lemma El except that we 
have to take care to choose S to be on the correct side of the contour g. Previously, it did not 
matter whether we formed S from the left or right components, and we arbitrarily chose the 
former. Now we choose S (either taking all the left or all the right components) in such a way 
that S H F>i }U = 0. This is possible because all the vertices in Bi )U are in a single connected 
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component (the contour g does not cross any edges whose endpoints lie in Bi v ). Now define 
a 9 as in the proof of Lemma [9] and continue as before. This establishes ([3]). 

For all 1 < s < r, and all u 6 U, there are at most 3 r x 4 contours V\ . . . v r for which u is 
on the left as we go from v s -i to v s : a choice of initial direction and direction at each step 
determines the contour. Summing over s and u, this implies that there are at most 4|C/|r3 r 
length-r contours with an edge whose dual connects vertices of U. By ([3]), 



p / cr u has a simple contour of length at least h which 
^contains an edge whose dual connects two vertices in U 



G v (Bi tV ) has parity-0 ones 



< 4\U\ r3 r max(/3A- 1/4 , 7 A 3/8 ) r 

r>h 

= 4|Z7|(3 x 0.238)' 1 J^(t + h)(3 x 0.238)* 
t>o 

. ., sh 3 x 0.238 + /i(l - 3 x 0.238) 

= 4|C7|(3 x 0.238) h — 



(1 - 3 x 0.238) 2 

238 

(1-3X0.238) 2 



There is a c > 1 such that 4(3 x 0.238)' 1 3x0 - 2 3 8 +M^ 3x 3 238 ^ < c '''for all sufficiently large h. 



□ 

By the upper boundary of C u we mean the set of all vertices of the form (x, v) for some x. 

Lemma 11. Let i G {0,1}. There is a c > 1 swc/i i/iai, /or all sufficiently large h and all 
v > h, the probability that a v has a simple contour that separates the set {—h + i, . . . ,h + i} x 
{0, . . . , h} from the upper boundary of C u , conditioned on the fact that cr u (Bi u ) has parity-0 
ones, is at most c~ h . 

Proof. Note that the separating contour cannot wrap around, owing to the boundary con- 
dition that a u (Bi^ u ) has parity-0 ones. If the separating contour has length r + 2 then its 
right-endpoint is in the set {(h + i + x — |, — |) | l<x<r}. There is a unique choice for the 
edge incident to each endpoint. Thus, there are at most r3 r possible contours. By Lemma [TUl 
the probability that cr u has such a simple contour, conditioned on the fact that cr u (Bi }U ) has 
parity-0 ones, is at most 



r&mxKtfX- 1 /*,^*)' 



\r+2 

r=h-2 

Thus, the probability is at most 



max^A-^^A 3 / 8 ) 2 £ r3 r max^A" 1 / 4 , 7 A 3 / 8 ) r , 

r=h-2 

which, as in the proof of Lemma [TUJ, is exponentially small in h. □ 

Lemma 12. There is a c > 1 such that, for all sufficiently large v , 

Pr(<r i/ has a cross contour ) < c~ u . 

Proof. Let g be a cross contour. There must be at least one other cross contour g' . For 
otherwise there would be a path p in C v from L(g) to R(g) such that <r u (V(p)) has parity-0 
ones or parity-1 ones, which would violate parity. Orient g and g' in opposite senses (one 
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away from V = —\ and one towards). Consider the connected components of the graph 
(y(C u ), E(C V ) \ {e* | e G g U g'}) , and let S be the union of all connected components that 
are left of either g or g' . Now proceed as in the proof of Lemma [9j using the fact that a cross 
contour has length at least u and the set of possible starting points has size 2u. □ 

Lemma 13. Let i 6 {0,1}. Fix v > 1. Conditioned on a v {Bi^) having parity-s ones (for 
any s € {0, 1} ), a u has no cross contour. 

Proof. A cross contour would have to cross a side edge in Bi v , which is impossible. □ 
Lemma 14. p = > . 

Proof. Fix v > 2. Suppose that (t v {Bq^) has parity-0 ones. If ov(0, 0) = then there is 
a simple contour of a u that separates (0, v) from (0,0). (Note that, by Lemma [TBI cross 
contours cannot separate these two vertices.) If the separating contour has length r + 2 then 
its right-endpoint is in the range (|, — 5), . . . , (r — 5,— 5). There is a unique choice for the 
edge incident to each endpoint. Thus, there are at most r3 r possible contours. By Lemma [TOT 
the probability that cr v has such a simple contour, conditioned on the fact that a u (Bo^ u ) has 
parity-0 ones, is at most 



^^maxOSA- 1 / 4 ,^ 3 / 8 )^- 

r=l 



Thus, 



Pr(ov(0,0) = I <t v (B Q)U ) has parity-0 ones) < max(/3A" 1/4 , 7 A 3/8 ) 2 ^ r3 r max(/3A- 1/4 , 7 A 3 / 8 ) r 

< (Q.238). 3 X 0.238^ 
~ v ; (1 - 3 x 0.238) 2 

< 1/2. 

Thus, p = > |. 

Similarly, suppose that a u (Bi jU ) has parity-0 ones. If there is no simple contour of o~ v 
that separates (1, v) from (1,0), then <r v (l, 0) = 0. We already saw that the probability that 
no such contour exists is greater than |. Thus, 

Pr(ov(l,0) = [ <T V {B\^) has parity-0 ones) > 1/2. 

So < |. Putting the two inequalities together, we have p^ < | < p = . □ 



4.7 In the absence of long contours, the spins of the terminals are nearly 
independent 

The ^-distance between two points (x, y) and (a/, y') in V(C U ) is max(|x — x'\, \y — y'\) where 
\x — x'\ is the minimum non- negative integer such that x = x' + \x — x'\ modulo 2v or 
x' = x + \x — x'\ modulo 2u. The ^-distance between two points (x, y) and (a/, y') in V(C*) is 
defined similarly. Let U x ^ be the set of vertices of C v whose *-distance from (x, 0) is at most h. 
Vertices in V(C U ) (or V(C*)) are *-adjacent if the *-distance between them is 1. A *-path 
on V(C U ) is a sequence v\, . . . , Vh of vertices in V{C V ) such that, for each j € {1, . . . , h — 1}, 
the vertices Vj and Vj+i are *-adjacent. A *-path on V(C*) is defined similarly. 
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Definition 15. Suppose v and h are positive integers. An /i-boundary of a vertex (x, 0) is a 
set of vertices B C V(C V ) such that the following are true. 

(i) (x,0) is not connected to (x,u) in the graph C u \ B, and 

(ii) B does not intersect U x ^u, and 
(Hi) B is a subset of U x h / 2 - 

(iv) The subgraph of C u induced by B is connected. 
See Figured 



(10,0) 




(0,0) 

Figure 5: Example /i-boundary of vertex (16, 0) with v = 20 and h = 5. 

Here is the relevant fact about /i-boundaries. If B is an /i-boundary of a vertex (x, 0) 
and cr(B) has parity-s ones then there is no contour of a which contains an edge whose dual 
connects two vertices in B. 

If B is an /i-boundary of vertex (x, 0) and B' is an /i'-boundary of (x, 0), then we say that 
B is inside of B' if every path in C„ from B' to (x,0) passes through B. Suppose that B 
and B' are /i-boundaries of (x,0) and, for some s G {0, 1}, o~(B) has parity-s ones and cr(B') 
has parity-(l ® s) ones. Then B n B' = 0, so exactly one of the following occurs. 

• B is inside of B' , or 

• B' is inside of B. 
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Definition 16. Suppose that s € {0, 1} and that k and d are positive integers. Let v = 2kd. 
Suppose that a is a configuration a : V(C U ) — > {0,1}. We say that a has phase s if the 
following holds for every terminal (x,0). 

• (x, 0) has a d-boundary B for which o~(B) has parity-s ones. 

• For every d-boundary B' of (x,0) for which cr(B') has parity-(KBs) ones, B' is inside 
ofB. 

Note that a configuration a can have exactly one phase (phase or phase 1) or it can 
have no phase. Suppose that v > 1 and that the configuration a : V(C U ) — > {0, 1} has 
phase s. Say that a d-boundary B of a vertex (x,0) is consistent if o~(B) has parity-s ones. 
Prom the consistent d-boundaries, we want to select a canonical one, that is in some precise 
sense "outermost" and also "minimal". For each terminal (x,0), let B x (o~) be the union of 
all d-boundaries of (x,0) which are consistent. Observe that B = B x (cr) satisfies the first 
three bullet points in Definition [151 but not the final one, as the subgraph C V [B] of C u 
induced by B may not be connected. Suppose that C U [B] has j connected components. 
Partition B x (a) = B x {a) U • • • U &{a) so that C u [B x (a)], . . . , C u [B x (a)] is an enumeration of 
these j connected components. Pach set B x (a) is itself a d-boundary. To see this, consider 
any vertex v E B x (a). From the construction of B x (cr), this vertex is contained in some d- 
boundary, which in turn is contained in B x (o~). So B x (a) satisfies the first three bullet points 
in Definition [151 i n addition to inducing a connected graph. 

From the first bullet point of Definition [15] it follows that B x (o~), . . . ,B x {o~) are nested; 
suppose the numbering indicates the level of nesting, with B x (o~) being the outermost. We 
now want to identify a minimal d-boundary within B x (o~). Let 

Ext Bl(a) = {v G V(C V ) | there is a *-path from v to (x, v) in V(C U ) \ Bl(a)}. 

denote the set of vertices lying in the "exterior" of B x (a). Finally define 

B x (a) = {v € B x (a) \ v is *-adjacent to some vertex in Ext B x (a)}. 

Note that B x (o~) is a d-boundary of (x,0) which is consistent. (Informally, B x (a) is the 
outermost such d-boundary.) To see this, observe that any path in the graph C v from (x,0) 
to (x, v) has a last vertex in the set B x (a), and this vertex must be in B x (a); this deals with 
the first bullet point in Definition [15l Consider the trail in the dual graph C* separating 
B x (a) and Ext B x (o~); the trail in the primal graph that shadows it at ^-distance \ inside 
takes in all the vertices of B x (o~) and establishes connectivity of C u [B x (a)]. This deals with 
the final bullet mark, and the remaining two are immediate. 

The d-boundary B x (a) is our desired canonical d-boundary, and it has the following im- 
portant property. If a' is a configuration that agrees with a on B x (a) U Ext £>;J.(<t), then 
Bx{o~') = B x (a). The reason is as follows. The set B = B x (a) is a d-boundary of (x,0) 
which is consistent with respect to a', i.e., cr'(B) has parity-s ones. It therefore gets incor- 
porated into B x (a') D B x (a) and hence into B l x (a') 2 B x (a). So Ext B x (a') = Ext B x (a) and 
Bx(o~') = B x {a). This fact becomes significant when we come to consider events supported 
on spins in the interior U = V{C U ) \(5U Ext-E?) of some d-boundary B. Specifically, condi- 
tioning on the event B x (a u ) = B (and on the phase s of er„) is equivalent to selecting cr u (U) 
according to the Gibbs distribution, with the boundary condition ll cr u (B) has parity-s ones". 
We refer to this property as canonicity of B x (a). 
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Before proceeding we need some definitions and observations concerning connected sub- 
graphs K of the graph (V(C*),a*). The *-diameter of K is the maximum, over pairs of 
vertices in K, of the *-distance between those vertices. We say that K reaches the lower 
boundary of C u if it contains a vertex of the form (x, — |) for some x. We say that it reaches 
the upper boundary of C v if it contains a vertex of the form (x, v + A). We say that it is a 
cross subgraph if it reaches both boundaries. We say that K wraps around if it contains the 
image of some path from (0, y) to (2i/, y) in Z x {0, . . . , u}, under the quotient map to V(C U ). 
We say that K is local if it is not a cross subgraph and does not wrap around. We say that 
K intersects U Xj h if some edge e* in K is dual to an edge e G E(C U ) with both endpoints 
in U X} h' A contour is a connected subgraph of (V(C*), a*) so all of the above definitions apply 
to contours. We refer to the connected components of (V(C*),o~*) as cr*-components. 

Lemma 17. If a has only local contours then a* has only local components. 

Proof. Suppose to the contrary that a* has non-local component K. Consider first the case 
of a cross component. 

So suppose K reaches both the upper and lower boundaries of C„, and that K has 2j > 
degree-1 vertices. (All vertices other than the degree-1 vertices have even degree, and the 
number of odd-degree vertices in a graph is even.) By adapting the standard algorithm 
for finding an Eulerian trail in a (connected) Eulerian graph, we may decompose K into 
j contours beginning and ending at degree-1 vertices. The method is as follows. Starting at 
a degree-1 vertex, trace out a trail in K subject only to the rule that we must turn through 
a right-angle at any degree-4 vertex. This trail can only end at another degree-1 vertex. The 
trail so formed is a contour; remove the trail from K and repeat j — 1 more times to obtain 
j contours in total. If any edges of K remain, start at any remaining vertex and trace out a 
closed trail that returns to the start vertex. Again, the rule is always to turn through a right- 
angle at any degree 4 vertex. Repeat until there are no edges remaining in K. We are left 
with j non-closed contours and an unspecified number of closed ones. Whenever a non-closed 
contour meets a closed one, we may splice the latter into the former, reducing the number of 
closed contours by one. Repeating as necessary, we obtain the sought-for decomposition of 
K into j contours beginning and ending at degree-1 vertices. 

If one of these j contours joins the upper and lower boundaries of C v we are done, as 
we have already found a cross contour and obtained a contradiction. Otherwise, there must 
be at least one vertex at which a lower-to-lower contour touches a upper-to-upper contour. 
Simply reroute the trails at this vertex to obtain two cross contours. 

Now consider the case where K wraps around. We may assume that K does not reach one 
of the boundaries of C v , say the upper one. Trace a closed trail along the upper boundary 
of K: this trail is a contour that wraps around, providing a contradiction. □ 

Since non-local contours are unlikely, Lemma [T7| allows us to concentrate on local a*- 
components. A o"*-component K that is local has a well defined inside and outside, and 
a boundary that is a valid contour. (More precisely, there is a canonical contour that has 
exactly the same edges as the boundary of K.) If K reaches neither the upper nor lower 
boundary of C„, then we may trace clockwise around K, always taking the leftmost option, 
until we return to our starting point. This procedure yields a simple contour; we refer to 
vertices of V(C U ) that lie within this contour as forming the interior of K, denoted Inti'T. 

If K reaches the lower boundary but not the top (or vice versa), then a slightly modified 
construction can be used. First lift K to a grid: for sufficiently large N, there is a connected 
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subset K of {1, . . . , N} x {0, . . . , v} which maps bijectively to K under the quotient map to 
C u . Note that lifting can only increase the diameter of K relative to K. We now have a 
natural ordering of the degree-1 vertices of K, namely by increasing x-coordinate. Start at 
the least degree-1 vertex in this ordering and and trace the boundary of K in a clockwise- 
leftmost fashion until the greatest degree-1 vertex is reached. This procedure yields a simple 
contour which partitions the vertices of V(C U ) into an inside and an outside (containing the 
point (0, y))\ again we refer to the former as the interior of K. 

Lemma 18. Consider a: V(C U ) — > {0, 1}. Let h < u be an integer multiple of 8. Suppose 

• a contains only local contours. 

• a contains no contour of length at least h/8 that intersects U x ,hl2- 

• a contains no simple contour separating U x ^u from the upper boundary ofC u . 

Then every a* -component that intersects U X; h/2 has *-diameter at most h/8. 

Proof. Suppose K is a <j*-component intersecting U x ^/ 2 . By Lemma [T71 K is local. The 
interior of K contains some vertex in U x ^i2- By assumption, the contour defined by the 
boundary of K does not separate U x ^ h u from the upper boundary of C u , so it does not 
separate U x ^/2 from the upper boundary of C v . The only remaining possibility is that this 
contour intersects U x ^/2, and hence has length at most h/8. It follows that K has *-diameter 
at most h/8. □ 

Lemma 19. Suppose that h is a sufficiently large multiple of 8, and v > h. Consider a 
configuration a: V(C U ) — > {0, 1} and a terminal (x,0). Suppose that the following are true. 

• a has no cross contours. 

• a has no simple contour of length at least h/8 that intersects U x ^/2- 

• a has no simple contour separating U x ^u from the upper boundary of C v . 
Then, for some s € {0, 1}, 

1. (x,0) has an h-boundary B for which o~(B) has parity-s ones. 

2. There is no h-boundary B' of (x,Q) for which cr(B') has parity-{\®s) ones. 

3. For any terminal (x',0) that has an h-boundary B' for which o~(B') has parity-s' ones, 
if a has no contour separating U x ^u from U x > then s' = s. 

4- If o~(Bo tU ) has parity-s' ones then s' = s. 

Proof. There are no contours that wrap around, since any such contour would either intersect 
Ugt fi/2, or would separate U x ^/4 from the upper boundary of C v . Thus, by Lemma [T71 all 
<7*-components are local. Let S be the set of all vertices in V(C U ) that are not in the interior 
of some <r*-component. That is 

S = V(C P ) \ U { Int K | K is a <7*-component}. 
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Note that a(S) has parity-s ones, for some s £ {0, 1}. Define S = V(C U ) \ S. Note that a(S ) 
in general has mixed parity; the salient feature is that cr(S) has consistent parity. 

We work first towards conclusion (1) of the lemma. By Lemma [TBI every o"*-component 
that intersects U x ^/2 has *-diameter at most h/8. Now (informally) we will construct the 
required /i-boundary by tracing round the inside of B xh / 2 , making a detour towards (x,0) 
around any cr*-components that stand in the way. (Recall that B x j l /2 = U x ^/2 \ U x u 2 -i i s 
the "goalpost" at distance h/2 from (x,0).) This strategy ensures we remain in the set S 
and, since all the cr*-components are small, our detours will not be too great. 

More formally, let W be the union of the set V(C U ) \ U x ^ h / 2 together with any sets of the 
form Int if that intersect B^. ^/2- That is, 

W = (V(C U ) \ U Xjh/2 ) U (J{Int if | if is a cr*-component and Int (if) n B xh/2 ^ 0}. 
The set 

dW = {v € V(C U ) | v $l W and v is *-adjacent to some vertex in W} 

is almost the /i-boundary B that we seek. Observe that any set of the form Int K is contained 
in a maximal set of the form Int K' , and the *-neighbours of Int K 1 are all in S. Thus dW is 
a subset of S, and necessarily has parity- s ones. 

And as we shall see presently, dW satisfies the first three conditions of an /i-boundary B 
- (i) every path from (x,0) to (x,u) intersects B, (ii) B n U Xj hu = 0, (hi) B C U x ^/2 ~ 
but not necessarily the final one, namely: (iv) the induced graph C^fS] is connected. (There 
may be islands of vertices of dW lying outside the /i-boundary we are trying to home in on.) 
However, we can ensure (iv) by defining B to be the subset of vertices in dW that can be 
reached from (x, 0) by a path in C u whose vertices all lie in V{C V ) \ W. 

For (i), observe that any path from (x,0) to (x, v) has a first vertex w in W. The vertex 
immediately preceding w is not in W but is adjacent to a vertex in W, and hence in B. 
(ii) follows from the fact that every vertex in Iffl U x ^/ 2 is within *-distance h/8 of a vertex 
in B x hj2- (iii) is immediate from the construction. To see (iv), denote by W° the set of 
all vertices in V{C U ) that can be reached from (x,0) by a path whose vertices all lie in 
V{C U ) \ W. Note that C U [W°] is connected and that B C W°. Let g* be the set of edges in 
C* separating W° and W; thus, e* £ g* iff e has one endpoint in W and the other in W°. 
Since C^fW 3 ] is connected, the edges in g* form a trail in C* starting and ending at vertices 
with y-coordinate — s. Following this trail anticlockwise, vertices in W° lie to the left and 
those in W to the right. In fact, the vertices immediately to the left of the £>*-trail (i.e., those 
at *-distance ^ from it) are precisely the vertices forming B: they are all *-adjacent to some 
vertex in W, and no other vertices in W° have this property. Thus, any two vertices in B are 
connected by a path, which is obtained by shadowing g* at *-distance ^. 

For conclusion (2) of the lemma, observe that there cannot be an /i-boundary B' of (x, 0) 
such that cr(B') has parity-(s©l) ones, as such a B' would have to exist entirely within S, 
and all *-connected components of S are small (*-diameter at most h/8). 

As for conclusion (3), it is impossible for s' ^ s. Consider the connected component 
of CjyfiS] containing B. If s' 7^ s then the boundary of this component contains a contour 
separating B and B', and hence U x hu and U^^u. If a(Bo jU ) has parity-s' ones for s' 7^ s 
then the boundary of the connected component of C^S*] containing B is a simple contour sep- 
arating B from the upper boundary of C v hence, separating U x ^u from the upper boundary 
of C„, establishing (4). □ 
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Corollary 20. Suppose that k > 1 that d is a sufficiently large multiple of 8 and that u = 2dk. 
Suppose a: V(C U ) {0, 1} has no contours of length at least d/8. Either a has phase or a 
has phase 1. 

Proof. Since there are no contours of length at least d/8 of any kind, the premises of Lemma[T9l 
are all satisfied for every terminal (x,0) and every other terminal (x',0). □ 

The following monotonicity property is useful for comparing different contours and bound- 
ary conditions. 

Lemma 21. Let x E TLj2vTL, let B be an h-boundary of (x,0) for some h, and let B' be an 
h' -boundary of (x,0) for some h! such that B is inside of B' . Then 

Pr(er„ has parity-0 ones at (x,0) \ ct v {B) has parity-0 ones) 

> Pr(<X;, has parity-0 ones at (x,0) \ cr u (B') has parity-0 ones). 

Proof. For each S C V(C U ) let a s : C v — > {0, 1} denote the configuration which has parity-0 
ones exactly on S. So o~ s (x, y) = 1 if and only if one of these two conditions hold: (x, y) £ S 
and x + y is even, or (x, y) ^ S and x + y is odd. 
For all X,Y C V(C V ) we have 

wc,M X )™cA<? Y ) = wcAa XnY )wcM XUY )(Pl) k 

where k is the number of edges uv G E(G) such that {(a x (u), o~ x (v)), (a Y (u), a Y (v ))} = 
{(0, 0), (1, 1)} (so either u € X and v ^ X and u £ Y and v € Y, or v £ X and u ^ X and 
v £Y and u € Y). In particular, 

^> x )^ c > y ) < wcM XnY )™cAo- XUY )- 

Let 

* = {S | {(x, 0)}UB'C5C f(C„)} 
y = {S\BlJB'cS<Z V(C U )} 

By the FKG inequality !)] we have 

(e^> 5 )) (e-c> 5 )) < f E ^(^)) f E 

Vse* / \sey J \seXAy / \sexvy J 

where X Ay is the family of sets X C FfC,,) such that ({(x, 0)} U B') n(5U B') = B' CX, 
and ^ V ^ is the family of sets X C FfC*,) such that {(x, 0)}UBUB'C X. Finally, 

Pr(cr^ has parity-0 ones at (x, 0) [ <r v {B') has parity-0 ones) = ^^ l5 ' e ' :y — ^i- ^ 

Pr(<r Jy has parity-0 ones at (x, 0) | cr u (B) has parity-0 ones) = ^^ e;t ' V ' y — ^ g . □ 

Lemma 22. There is a c > 1 suc/i that the following is true for any k > 1, any s € {0, 1}, 
any sufficiently large d which is a multiple of 16, and any assignment {B x } of d-boundaries 
for each terminal (x, 0): 
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• For every parity-s terminal (x,0), 

| Pr(<x M (x, 0) = 1 | a k ,d has phase s and B x (cr kjd ) = B x ) - p = \ < c~ d . 

• For every parity-(l@s) terminal (x,0), 

| Pr((r M (x, 0) = 1 | <x M has phase s and B x (a k4 ) = B x ) - p^\ < c~ d . 

Proof. By symmetry (rotating the gadget so that parity-0 vertices become parity-1 vertices 
and vice- versa), it suffices to prove the inequalities for s = 0. For any m > 1, define 

p~(m) = Pr(er m (0, 0) = 1 | <x m (.Bo,m) nas parity-0 ones), and 
p^(m) = Pr(er m (l,0) = 1 [ cr m (Bi tm ) has parity-0 ones). 

Now note that for any v > m, p = (m) and p^(m) (as defined above) are the same as the 
equivalent expressions in the gadget C u . In particular, 

p = (m) = Pr(<j Jy (0, 0) = 1 | <X;,(i?o,m) has parity-0 ones), and 
p^(ra) = Pr(<j Jy (l,0) = 1 | cr u (Bi :m ) has parity-0 ones). 

Thus, by fixing large v and increasing m, Lemma [21] implies that p = (m) is weakly 
decreasing in m and that p^(m) is weakly increasing. Thus, p = = lim m _>. 00 p = (m) and 
p^ = lim m _ i , 00 p^(m). Also, for a parity-0 terminal (x,0), the target probability 

Pr(cr M (x,0) = 1 | er k ,d has phase and B x (cr ktd ) = B x ) 

is between p = (d/2) and p = (d / 4) . Similarly, for a parity-1 terminal (x, 0), the target probability 

Pr(o- fejrf (x, 0) = 1 [ a k ,d has phase and B x ((T k ,d) = B x ) 

is between p^(d/A) and p^(d/2). (Here we use crucially the canonicity property of B x { ); refer 
to the discussion following Definition [161) Thus it suffices to show 

P = {d/4:) <p = + c~ d and p^{d/4) > p* - c~ d , 

First we take a qualitative step. Pick w > 8d sufficiently large that p = (w) < p = + 
I ^o,d/4 1 (c') — rf / 16 and p^(w) > — \U ltd / 4: \(c')~ d ^ 16 , where d is the maximum of the constants 
given in Lemma [10] and Lemma [11] This can be done since d is sufficiently large and p = = 
lim m _>oop = (m), and p^ = linim^oo p^ (m) though w may be quite a lot larger than d. 

For i € {0, 1}, let Fj be the event that there is a <i/2-boundary B of vertex (i, 0) in gadget 
C w such that cr w (B) has parity-0 ones. Recall from the definition that a (i/2-boundary of (i, 0) 
is a subset of U^u. Let Ei be the event that a w (Bi^ w ) has parity-0 ones. 

For each i € {0, 1}, applying Lemma [131 and Lemma [lOl with h = d/16 and U = U^ d u, 
and Lemma QT] with h = d/8, we find that, the conditional probability that the following 
hold, conditioned on Et, is at least 1 — 2\U i ^ d / /l \{c')~ d / l& . 

• <r w has no cross contour. 
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• cr w has no simple contour of length at least d/16 which contains an edge between two 
vertices in U^u. 

• cr w has no simple contour that separates Ui^/s from the upper boundary of C v . 

Now, applying Lemma [19] with h = d/2 and v = w and x = i, if all of these hold and event 
Ei occurs then event occurs. Thus, 

Pr(^|^)>l-2|C/ M/4 |(c')- d/16 . 

But by Lemma EH we have 

Pr(<x„,(0,0) = 1 | F AE ) > p = (d/4), and 
Pr(<r„,(l,0) = 1 | F 1 AE 1 )<p*{d/A). 



So 



p=( d / 4 ) < Pr(cr w (0,Q) = IAF \E ) < Pr(<r tt (0,0) = 1 \ E ) _ p=(w) 



Pr(F | Eo) - Pr(F | E ) Pr(F | E ) 

- 1 -2i^lr^ -^ ( '" )+4|Cfwl(c ' rd/ ' 6 - 

since 2 1 1/ ^ /4 1 (c')~ rf/16 < 1/2. 

A similar inequality holds for p^(d/4): 

, Pr( g|l ,(l,0) = lAF 1 |g 1 ) 

^ 1 7 ' ~ Pr(Fi | £7i) 

> Pr(cr„(l, 0) = 1 | E{) - Pr(^i | ^) 

>^M-2|c/ M/4 |( c ')- d/16 . 

Thus, 

P=(d/4) <P = H + 4|C/ 0j(i/4 |( C ')- d/16 <P= + 5\U 0>d/ ,\(cT d/W , and 

p^w) > p^h - 2\u hd/4 \(c')- d/w >v*- mi, d /4\(c'r d/w - 

The result follows by noting that [E^mI is 0(d 2 ) and picking c = (c') 1 / 17 , say. □ 

We now prove the main proposition. 

Proposition [5], There is a c > 1 such that, if d is a sufficiently large multiple of 16, k is an 
integer greater than or equal to 1 and r is a configuration r : d — > {0, 1}, then 

I Pr(<x M (T M ) = r) - Mm (t)| < c^fe 2 . 

Proof Fix fc > 1, d a sufficiently large multiple of 16, and r : — > {0, 1}. Let c' be the 
minimum value of the constant c from the Lemmas [UJ [T21 and [221 

The probability P^o-fc^T^) = t) is the sum of the following probabilities (conditioned 
on disjoint events) 

• Pr(<T/ c c ;(Tfc j rf) = t | ak,d does not have a phase) Pr(cr fc rf does not have a phase). 
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• (summed over all assignments B x of d-boundaries for each terminal (x,0)) 

P^(ck,d(Tk,d) = t I (Tk,d has phase and for all terminals (x, 0), B x ((Tk,d) = B x ) x 
Pr( (Tk ; d has phase and for all terminals (x,0), B x (ak,d) = B x ) 

• (summed over all assignments B x of (i-boundaries for each terminal (x, 0)) 

P r (°"A;,d(TA;,d) = 7" I < J k,d has phase 1 and for all terminals (x,0), B a i( 'fe ) d) = ^s) x 
Pr( <jfc^ has phase 1 and for all terminals (x, 0), B x {a^,j) = B x ) 

By Lemmas [9] and [12] and Corollary [20] the probability of the first of these is at most 
2\V{C* k4 )\{c')- d / % . (We will use this below.) 

Now consider an assignment B x of (i-boundaries for each terminal (x, 0). For any two 
terminals (a/,0) and (x",0), the random variables <Xfc )£ j(a/, 0) are tTk,d(x", 0) are independent, 
conditioned on the fact that (Jk,d has a given phase, and for all terminals (x, 0), Bx^tr^d) = B x . 
Also, by Lemma l22l for all s € {0, 1}, 

• For every parity-s terminal (a/,0), 

|Pr(<Tfc rf(x', 0) = 1 | cr/j d has phase s and for all terminals (x,0), B x {cr k4 ) = B x )—p = \ < 

• For every parity-(l©s) terminal (a/,0), 

| Pr(crfc i( i(x / , 0) = 1 | (Tk,d has phase s and for all terminals (x,0), B x (crk,d) = Ac) - < 

Now, for any probabilities a\, b\, . . . , a^, bk, we have 

k k 



i=l i=l 



k 

Ql ■ ■ ■ aj-i( a j - bj)b j+1 ...b k 



i=l 



so if we fix a given phase, and r assigns spin 1 to A/ terminals whose parity agrees with that 
phase, and spin 1 to k" terminals whose parity disagrees with that phase then, letting 

P = (P = ) fc '(! ~ P = ) k ~ k> {P^) k (1-J^) fe k , 

we have 
p-2k{cT d 

< Pr(cr/ c ^(Tfc ]( i) = r | cr/c^ has the given phase and for all terminals (x,0), B x (crk,d) = B x ) 
<p + 2k(c')- d . 

So summing up, Pr(<r fe i( i(Tk,d) = t I &k,d has the given phase) is between p — 2k(c')~ d 
and p + 2k(c')~ d so, since the phases are equally likely, 

Mfc,d( r ) ~ 2k(d)~ d < Pr(<x M (T M ) = r | <r M has a phase) < Mm( t ) + 2k(c')~ d . 

Finally, since the probability that <r fcjd has no phase is at most 2\V(C k d )\(c')~ d ^ 8 , as we 
observed above, 

M M (r) - 2(* + |y(^ d )|)(c')^ /8 < Pr(<r M (T M ) =t)< /* M (r) + 2(fc + |^(C^)|)(c')- d/8 . 
The proposition follows by choosing c to be sufficiently small with respect to d . □ 
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5 Proof of Theorem [T] 

5.1 Efficiently approximable reals 

Lemma 23. Suppose that f3, 7 and A are efficiently approximable reals satisfying ([1]). Then 
p = and are efficiently approximable reals. 

Proof. Recall that p^ > (LemmaH]). Let q be a multiple of 16 greater than (2+log 2 (l/p^))/ log(c) 
where c is the constant given by Lemma [221 Consider the following algorithm. 

• Input an error parameter < e < 1/2. 

• Set m = q\log(e^ 1 )~\ . 

• Compute rational approximations /3, 7, A satisfying 

p e -e/8\E(C m )\ < p < ^ e e/8|S(C m )| 
-s/8|B(C m )| < ~ < 7e£ /8|i?(C m )| 
Ae - £ /8|y(C m )| < A < Ae e/8|K(C m )|_ 

• Using the algorithm of [27\ Theorem 2.8], compute 

Z = ^pb(<r)z ( c{a)yi{<r) 

(7 

Z> = Y fi b< y°) 
o-:cr(0,0)=l 

ct:<t(1,0)=1 

where the sums range over configurations <r of C m such that cr(-E>o,m) has parity-0 ones. 

• Output Z'/Z as the approximation to p = , and Z"/Z as the approximation to p^. 

For the computation of Z, Z' , and Z" we use the fact that the grid graph C m \ i?o,m has 
treewidth m [H Corollary 89]. We also use the fact that its tree decomposition is easy to 
compute. So this algorithm runs in time bounded by a polynomial in 1/e. We will show that 
the algorithm is an FPRAS for p = and p^. Define 

w = ^2^)^)^) 

a 

w / = J2 p^)^)^) 

cr:o-(Q,0)=l 

W" = /3 6(<t) 7 c(ct) A £(ct) , 

a:a{l,0)=l 

where the sums range over configurations a of C m such that cr(i?o,m) has parity-0 ones. 
For any a we have 
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This implies e £ l i W < Z < e £ '^W and similarly for Z' and Z", and therefore e £ I 2 W' /W < 
Z'/Z < e £ / 2 W'/W and e- £ / 2 W"/W < Z"/Z < e £ / 2 W"/W. We will show 

p= < W'/W < p=e £/2 (4) 
e" £/ V < W"/W < p + (5) 

W'/W and W" /W are just the probabilities that an even or odd terminal gets assigned 1 
in a random configuration of C m , conditioned on a certain 2m-boundary. By Lemma [2T1 we 
have p= < W'/W and W" /W < p+ for any m, establishing the first inequality in ([3D and the 
second inequality in (J3J). 

By Lemma [22j there exists c > 1 such that 

W'/W <p= + c - qlos ^ =p=(l + e qlo ^ c) /p=). 

Since 

£ ( 9 log(c)-l) < /]/2N(glog(c)-l) < 
which is less than p = by Lemma [U we have 

W'/W <p = (l +e) < e £ . 
This establishes ([!]). Similarly, by Lemma [22l 

w"/w >p + - c - qlo ^ V) 

> p^{l-e qlo ^ c) /p*) 
>p*{l-e/2) 

> p^e~ £ 

This establishes (J3J). □ 

Lemma [23] gives us a way to obtain multiplicative approximations p = and p^ of the real 
numbers p = and p^ . When we use these approximations, we will need to know that 1 — p = 
and 1 — p^ are also good multiplicative approximations to 1 — p = and 1 — p^, respectively. 
As we show below, this follows from the fact that p = and are in (0, 1) (which follows from 
Lemma U and Lemma I14|) . The following lemma gives us what we need. The reason for 
introducing the rational p' in the statement of the lemma is that, since it is rational, it can 
be hard- wired into any algorithms (whereas a real number can't be). 

Lemma 24. Suppose that p G (0, 1) is an efficiently approximable real number. Let p' be 
a positive rational with p < p' < 1. For any 5 £ (0,1), and any real number p satisfying 

e -5(l-p>)/2~ < p < e «(l-p0/2£ we have e -6(i -p) <\-p< e S (l - p). 

Proof. Let 5' = 5(1 — p')/2. Since p > e~ s 'p > p(l — 5') > p — 5' and similarly p >p — 5' , we 
have 

P-ri(i-|^)-(i-rt-^i-^0-rt^-0-rt(i + |^ 

Thus, 

(1 - p)(l - 5/2) < 1 -p < (1 - p)(l + 5/2), 
which suffices. □ 
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The following problem is NP-complete |12| . 



Name PlanarCubicIS. 

Instance A planar cubic graph G and a positive integer h. 

Output "Yes" , if G contains an independent set of size h, and "No" , otherwise. 

Suppose that (3, 7 and A are efficiently approximable reals satisfying ([1]). We will give a ran- 
domised polynomial-time algorithm for PlanarCubicIS, using as an oracle, an FPRAS for 
DegreeFourPlanarTwoSpin(/3, 7, A). The oracle will be used to approximate Z l -~ X (G), 

for some suitably-defined parameters 7 and A, where 7 is exponentially small in and A 

is exponentially large. From this, it will be easy to determine whether G has an independent 
set of size h. 

Lemma 25. Suppose that f3, 7 and A are efficiently approximable reals satisfying ([1]). There 
is a polynomial-time randomised algorithm that, given a planar cubic graph G with \V(G)\ 
sufficiently large, outputs planar graphs J and J' with maximum degree at most 4 and ran- 
domised approximation schemes for positive reals K , 7 and A. The running time of each of 
these approximation schemes is bounded from above by a polynomial in \ V{G)\ and the desired 
accuracy parameter e. With probability at least 14/15, the parameters satisfy A > 4^( G ) and 
7 < \-\ y ^\ and 

e-y%^{G) < < e^Z^G). (6) 

Proof. Let G = (V,E) be a planar cubic graph and let n denote \V\. 

The algorithm for constructing J and J' uses a quantity 5 € (0, 1). It will be important 
for the proof that 5 is sufficiently small. Rather than giving a technical definition here, we 
introduce upper bounds on 5 in natural places throughout the proof. The reader can verify 
that each of these upper bounds is at least the inverse of a polynomial in n (so the algorithm 
runs in polynomial time). 

The first step is to use the given FPRASes for /3, 7 and A, and the FPRASes for p = and 
p^ from Lemma [23] to compute values j3, 7, A, p = and p^ satisfying 

e" 5/3 /3 < $ < e 5/3 (3, 
e -<V3 7 < ^ < e <5/3 7j 

e- 5 / 3 X < A < e 5 / 3 A, 

> = , 

(7) 



e- 5 /V<P^<e 5/ V 
e- s/3 (l-p=) < 1 -p < e s/3 {l-p=), 
e -«5/3(! _ p ?) < 1 _ p < e «/3(! _ 

P>1. 

The first five lines in ([7]) follow directly from the definition of FPRAS in Section [3l The next 
two lines follow from Lemma [24"1 using the fact that p = and p^ are in (0, 1), as argued just 
before Lemma l24l Since j3 > 1 by (pQ), we can ensure that (3 > 1 by taking /3 to be the 
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maximum of 1 and the output of the FPRAS. For this step we adjust the failure probability 
of the FPRASes (as described in Section [3|) so that the probability that Equation ([7]) fails to 
hold is at most 1/15. Note that the running time of the FPRASes is polynomial in 1/5 (even 
though the application of Lemma [23] means that we have to call the FPRASes for p = and 
with slightly smaller values 5' .). 

We will show below how to use G and these approximations to define positive integers 
k\, k>2 and d, which will be used in the construction of J and J'. These quantities will be 
bounded from above by a polynomial in n. 

We first show how to construct J and J 7 , using k%, A>2, d and k = max(A;i, 3A?2)- The 
high-level construction is illustrated in Figure El 




Figure 6: An illustration of how G is transformed into the graphs J and J' . A fragment of G 
is shown on the left. The graph J is a collection of copies of C k} d, one for each vertex of G. 
A fragment of J is shown in the middle. The copies of C k d are shown as grey annuli. The 
corresponding fragment of J' is shown on the right. The stripes represent the sets of edges 
between copies of C k ,d m J'- J' also contains some "bristles" (described later) which are not 
shown. 

The construction of J is straightforward. Essentially, J consists of |V| copies of Ck, d , 
with one copy for every vertex in V . Thus, the vertex set V(J) is the set of ordered pairs 
V(J) = V x V(C k)d ) and the edge set E(J) is given by E{J) = V x E{C k4 ). We will use 
C[u] to denote the gadget corresponding to vertex u G V. Formally, C[u] is the graph with 
vertex set {u} x V(Ckd) an d edge set {u} x E{Gj. ( j). To simplify the notation, for u € V 
and < j < k — 1, let T l [u,j\ denote the j'th parity-1 terminal of C[u\. Formally, this is the 
vertex (u, (Ajd + 1,0)) of J. Similarly, let T°[u,j] denote the j'th parity-0 terminal of C[u\. 
Formally, this is the vertex (u, {^jd + 2j, 0)) of J. Let T[u] be the set of terminals of C[u]. 
Let flu and \i u be the distributions on configurations a : T[u] {0, 1} corresponding to 
the distributions d , n\ d and \i defined in Section |4~H 

To simplify the description of J', consider a planar embedding of G in which each vertex 
u of G is associated with three "endpoints" no, u\ and U2, which are arranged together in 
clockwise order in the plane. The edge set E can then be viewed as a matching M on the 
points LLgyi^O' Ul > u<2 ^ sucn that 

• (u,v) E E if and only if there are exactly two points Ui and Vj such that (u{,Vj) G M, 

• No two edges of M. cross. 

The vertex set V(J') consists of V ( J) , together with a set of nk\ new vertices, called "bristles" . 
Formally, V(J') = V(J) U{(u, j) \ u G V, < j < k± — 1}. Finally, the edge set of J' consists of 
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E(J), together with new edges connecting the bristles to the parity-1 terminals of the gadgets, 
and new edges matching the parity-0 terminals of the gadgets (guided by the matching M). 
The edges connecting the bristles to parity-1 terminals of the gadgets are those in the set 

E B = {{{u,j),T l [u,j}),u G V,0 < j < h - 1}. 

It is more complicated to describe the edges matching the parity-0 terminals of the gadgets. 
The idea (see Figure[7|) that if u a is matched to v b in M. (where a G {0, 1, 2} and b G {0, 1, 2}) 




Figure 7: Terminals of u are matched to terminals in v, reversing the order. 

then the parity-0 terminals T [u, ak 2 ) , . . . , T° [it, ak 2 + k 2 — 1] get matched to the parity-0 
terminals T°[u, bk 2 ], ■ ■ ■ , T°[u, bk 2 + ^2 — 1]- However, there is a further complication: To 
ensure that J' is planar we must ensure that one of these sequences of terminals is matched 
in clockwise order, and the other in anti-clockwise order. Thus, let 

Em = {{T°[u, ak 2 + j], T°[v, bk 2 + k 2 -l- j}) \ u<v, (u a ,v b ) G M, < j < k 2 - 1}. 

Then E(J') = E{J) U Eb U Ej^. Note that both J and J' are planar as required. 
We next show how to define the positive integers fei, k 2 and d. Define 



P 



1 — p p 



M P \ 5 X M P' 



W = P 



P 1 
1 7 



where P t denotes the transpose of the matrix P. Also, define 



P 



I — p- p- 
l-pt p^ j 



.1 [ - P [ ^ l_) P. 



W = P 



P 1 
1 7 



Note that if © holds then, for any s G {0, 1} and s' G {0, 1}, 



e S P S , S ' < Ps,s' < eS Ps,s> 
e~ 5 M S:S , < M S)S , < e 5 M s 

e' 5 W s <W S < e 5 W s . 



(8) 
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The matrix M has the following informal interpretation. Suppose that two parity-0 ter- 
minals t and t' are adjacent in J', and that a: V(J') — > {0,1} is a configuration. If these 
two terminals have spins a(t) and a(t'), respectively, then the edge between them contributes 

a factor ( f i ) to Wj>{cr). We will show below that, if t is a terminal of C[u] and 

the spins of C[u] are chosen from the idealised distribution p^, then the probability that the 
spin of terminal t is j is P s j- Thus, informally, M s s i captures the expected contribution of 
this connection (in the idealised distribution), where s' represents the phase of the gadget of 
terminal t' . 

The informal interpretation of W is that, given any configuration a : V{J) —> {0, 1}, a 
parity- 1 terminal t which is connected to a bristle b will contribute a factor to 

the sum ]P , wji(cr'), where the sum is over all configurations a' : V(J') —> {0, 1} which agree 
with a except possibly at the bristle b. This informal description is just to provide intuition 
— the technical details are given below. The main idea is that, if the spins of the terminals of 
the gadgets are chosen from the "idealised" distribution then, if the gadget of t has phase s, 
then the terminal t will contribute a factor of Wi® s to the expected contribution from this 
bristle. 

We now introduce some calculation which will be needed to describe the algorithm's 
computation of k%, A>2, and d and also to give the definitions of the real numbers 7 and A. 
The first step is deriving some tedious but necessary bounds on the various quantities defined 
above. In particular, we will define positive rational numbers A - and A + and a rational 
number £ 6 (0,1) (independent of 5 and n, but depending on /3, 7 and A). These will be 
hard-wired into the algorithm. We will prove that, provided that 5 is sufficiently small, each 
M Sj8 i and W s satisfies A" < M SyS > < A+ and A~ < W s < A+. Also, each of - p?, 1 - 7, 

p£ and A is at least £. We will also prove that p£ < 1 (and, from ([7]), we have (3 > 1.) Finally, 
we prove [3rj < 1 — £. Here are the details (which the reader may skip). 

• By Lemmas H] and Q31 we can define positive rational numbers p~ and p + such that 
every element P SjS ' of matrix P satisfies p~ < P s s i < p + . Then 

{p- f{2 + (3 + 7) < M s , s , < {p + f{2 + p + 7), 

so, since 5 < 1, 

e- 1 (p-f(2 + (3 + 7) < M SiS , < e(p + f(2 + /3 + 7), 

so to get the required bounds, we take any A~ < e~ 1 (p~) 2 (2 + j3 + 7) and any A + > 
e(p+) (2 + p + 7). The bounds on W s are similar. 

• To ensure p = — p£ > £, choose rational numbers pi, p2, p3, and ^4 such that p^ < p\ < 
P2 < P3 < Pa < P = - These exist by Lemma [TH which guarantees that p^ < p = . Then, if 
5 < P4 — P3, Equation (|7|) guarantees p = >p = e~ s >p = {l — 5) > p = — S > p = — (/04 — P3) > 
ps. (Note that the calculation used p = < 1.) Similarly, if 5 < (p2 — Pi)/2, Equation ([7]) 
guarantees p£ < e s p^ < p^(l + 25) < p^ + 25 < p^ + (p2 — pi) < P2- (Again, we used 
p^ < 1.) It suffices to take any £ < p% — P2- 

• We can similarly establish 1 — 7 > £ by considering a sequence of rational numbers 
between 7 and 1 (using the fact that 7 < 1 by (p}) and we can establish p^ > £ by 
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considering a sequence of rational numbers between and p^ (using the fact that p^ > 
by Lemma HI). 

• Then, by JT]), A > 0, so taking any £ < A, we can choose a rational number £' with 
<£<£'< A. Then choosing 5 < log(£'/£) ensures e -5 ^ > £ so A > e^A > e _<5 £' > £. 

• Note that the second bullet point already establishes p^ < 1. 

• Finally, ([1]) guarantees /?7 < 1, so choose rationals (3' > (3 and 7' > 7 with /3'7' < 1. 
Choose £ sufficiently small that j3'y' < e -3 ^. Then choose 5 < £/2 to ensure 

£7 < e 2,5 /37 < < e~ 2 « < 1 - f . 

We can make the following conclusions. 



Mi,! - M ,i 


= (p 


-p^)((/3- 


i)(i-^) + (i-tV) >e 3 


M 0) i - M ,o 


= (p 




i)(i-P) + (i-7)P) >e 3 


wi-wo 


= (p 




i) + (i-7)A)>e 3 - 



We can now define &2. Since 

M ,oMi,i - Mq^ = det(M) = det(P) 2 (^7 - 1) = (p - ?) 2 (^y - 1) < ~f, so 



Mq^Mx,! = Mo^ + CMo.oMLi-Mo^) < M 2 ^-g 3 = 1 _ < x _ _J^_ < g - € 3 /(A+)2 



^0,1 
Then let 



^0,1 



Mq X ~ (A 4 



A: 2 



(ra 2 + n)21og(5)(A+) z 



Then, if we ensure that 5 < (£ 3 /(A + ) 2 )/8, we have 

k 2 



4S M 0fi M hl 



< e 



-k 2 e/(2(A+) 2 ) 



< 5" 



Then define 



By 



7 



k 2 



< e , Sk2 / Mo,oM ltl 

V m 2 . 



< 5" 



Also, there is a randomised approximation scheme for 7 whose running time is at most a 
polynomial in n and in the desired accuracy parameter e. 
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Next, we will define k\. Recall that W\ > Wq and M\\ > Mq^. If n is sufficiently large 
with respect to log(A + /A~) > log(Wi/Wo) then there is a positive integer k\ (which the 
algorithm can compute) satisfying 

3A&tog(M u /M ,i) , log(4.1)n 3fc tog(Mi,i/M ,i) , Iog(4.9)n 
H — — < «i < — — h 



logiWx/Wo) log(Wi/W ) log(Wi/W ) logCWi/Wb) 

Note that h = 0(n 2 ). Also, 



(4.1) n < 

and 



'- — - \ Jfcl / - — - \ 3k2 

Wo) I Mi J 



Wo) VMiJ 



< (4.9) n . 



Now if we ensure 5 < n 2,5 then, for sufficiently large n, 5 < n log (4. 1/4)/ (2 k i + 6/02) and 
5 < relog(5/4.9)/(2fci + 6fc 2 ), so 

4 . < (V-Ef (.-SO"" 

V Wo I V 



and 



Note that 



e 



- — • \ fei / - — - \ 3^2 



5 n < ' 



fcl / i» \ 3fc2 



Then define 

s_r. , , 

Note that there is a randomised approximation scheme for A whose running time is bounded 
from above by a polynomial in n and the desired accuracy parameter e. Also, A < so 

7 < A~ n and A > 4 n , as required. 

Now let k = max(fci, 3&2)- Finally, the gadget will use a parameter d. By Proposition El 
there is a c > 1 (not depending on k) such that, for all sufficiently large d which are multiples 
of 16, and all configurations r : T^^d — > {0, 1}, 

I Pr(<x M (T M ) = r) - Mm (t)| < c^fc 2 . 

The algorithm will choose d to be a multiple of 16 such that d = 0(n 3 ) and 

jc- — - , k\n , r- — • .k2\E\ 

c -d k 2 < (e~ S W ) (e- s M 1A ) 
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This can be done, since \E\ = 0(n). We will use below the fact that 

w km M k2\E\ 

max I Pr(<ji, ATu d) = r) — iiu At~)\ < ; ■ , (9) 

r:T fe , d ->{o,i} 1 1 ' y ' ' ' " ' V n (/3 + A) fcin /3 fc ^l2 2fc ™2™n 

which follows from Equations (|7|) and ([8]). 
Let K be the positive real given by 



K 



2" 



W kin M k2\E\ 



Note that there is a randomised approximation scheme for K whose running time is at most 
a polynomial in n and in the desired accuracy parameter, e. 

All that remains is to establish ([6]), which we do in the remainder of the proof. Let 
T = U U £yT[u] be the set of terminals in V(J). For every configuration r : T — > {0,1}, let 
wt(r) = Yla<=V(J)^{o i}-a(T)=T w j( a )- The quantity wt(r) is the contribution to Zp a> x(J) from 
configurations a with cr(T) = t. Similarly, let wt'(r) be the contribution to Zp^^\(J') from 
these configurations. Since V(J') = V(J), we have wt'(r) = Effev(j)-).{o i}'<j(t)=t ( ff )- 
Let F(t) denote wt'(r)/wt(r). Then 

Z(*,yM') = £^ { o,i } wt'(T) = E^ { o,i } wt(r)F(r) = E[F(<jj(m (1Q) 
We can write F(t) in terms of j3, 7, and A: 



F(r) 



nn 



u&V j=Q Vvl 7 ^ V A // r (T>,j]) 



k 2 -l 

n n 



$ 1 

1 7 



(u a ,v b )eM,u<v j=0 V f/ T(T°[u,ak 2 +j]),T(T°{v,bk 2 +k 2 ~l-j]) 

(11) 



We now define an "idealised" distribution on configurations assigning spins to the termi- 
nals. First, let <x be a random variable which is drawn uniformly from V — > {0, 1}. Each 
realisation a of a can be thought of as specifying, for every vertex u € V, a phase cr(u) G {0, 1} 
for the gadget C[u]. Conditioned on the realisation a = a, the random variable <x : T —> {0, 1} 
is distributed as follows: for each u € V, cr{T[u\) is chosen independently from the distribu- 
tion im}^ . Note that the (unconditioned) random variable cr has the property that <r(T[?x]) 
is distributed as independently of all <r(T[u']) for u' 7^ u. 

We wish to estimate E[i ? (crj(T))], but the distribution of aj{T) is somewhat complicated. 
Instead, we will first estimate E[i ? (<r)], and we will later use Proposition [5] to show that these 
two quantities are close. From the definition of <r, we have 

"W)] = ^ E E[F(<x)|<r = a]. 

5-:V-!-{0,l} 

Now we will argue that if <j(u) = s then, for every parity-0 terminal t of C[u], it is the 
case that Pr(<r(t) = j \ a = a) = P s j. (To see this, consider the possible cases. If s = 
then, from the definition of ^ d> the probability that &(t) = 1 is p = , which is Po,i) but the 
probability that &(t) = is 1 — p = = Pq,o- The situation is similar if s = 1.) On the other 
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hand, similar reasoning shows that, for every parity-1 terminal t of C[u], it is the case that 
Pr(<x(i) = j | cr = a) = P\® B j. Thus, we have 

1 - 

E [ F (°)} = ^ E Y F ( T )YlYl P a(u),T(TO[u,j])Pim(u),r(T^[u,j])- 

ct:V->{0,1} t:T->{0,1} M€Vj = l 

Plugging in (fTT|) . the contribution of each a to the right-hand-side of the above equality 
is 2~ n multiplied by the product of the following terms: 

nn e ' 

«gy j=o se{o,i} v v 7 v 

fe-i 

n n e p ie^w^' 

uev j=k! se{o,i} 

n n e e Pa(u), s Pa(u),s' L i) , 

M££ i=o se{o,i}s'e{o,i} v '/«,«' 

fe-1 

n n e 

uGVi=3fc 2 se{o,i} 



The second and fourth of these terms are equal to 1, and the first and third simplify using 
the matrices that we defined earlier, so we get 



nF(^)} = ^ Y U w i%hu) II Mtwm 

ct:V->{0,1} u&V (u,v)&E 

Then, plugging in our notation from earlier, we have 



kz 



2 n 



E 

:V->{0,1} 



Mo, 



Replacing cr(it) with 1 © &(u), we get 

J:V^{0,1} 



Mi.i 
M ,i 



/ Mq,q 

VM i 



fc2c(<r) 



my 
w y 



f-1 



fei£(cr) 



Since G is cubic, we can count the pairs (v, (u, v)) £ V x E with = 1 in two ways to 
get U{a) = 2c{a) + (\E\ - b(a) - c(a)). So b{d) = c(a) + \E\ - 3£(a), which implies that 



w£ in M* 2 l El 



Mo 7 



E 

cr:V-i-{0,l} 



M 2 i 



k 2 



WA kl fM ^ 3k2 



W J VMi 



(12) 
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Plugging (fTU|) and (fT2"j) into ©, it remains to prove 

e" 1/4 E[F(6-)] < E[F(aj(T))] < e 1/4 E[F{&)}. (13) 
Let tp = \E[F(aj(T))] -E[F(&)]\. Now 

t/> < ( max F(r)] £ |Pr(^j(T) = r) - Pr(<r = r)| . 

To emphasise that summation over t : T — > {0, 1} can be broken into summation over each 
restriction t(T[u]), we will write the summation index as Vu, t(T[u\) : T[u] — > {0, 1}. By (jTTj) . 
^(r) < (/3 + \) kin ^\ E \ 7 so we can write 



V-u,r(T[u]):TH->{0,l} 

= (/3 + A) fein /3* 2|B| 

Vit,r(T[w]):TM-*-{0,l} 



J] Pv(aj(T[u]) = r(T[u])) - ]J Pr(<x(T[u]) = r(T[u])) 
ugv uev 

J] Pr(*rj(T[«]) = r(T[u])) - [] M(r(T[u])) 



Using the inequality | Y\ «« — n &«l < S l°« ~~ M va hd for values < a u , 6 U < 1, as in the 
proof of Proposition [5j we have 



V ,< (/ 3 + A )^ 2 |^l 2 ^ n ma x |Pr(<7 M = T)- MM (r)|. 



Applying 



^ < (r + X ) k ^f3^ E \2 2kn n =2 hi = _J y_ = 

V ^ I (/3 + A) fcl "/3 fc 2l s l2 2te 2™n; 2" 



So to establish (fT3|) , we note that 

1 1 ELF(<xj(T))l 1 1 

1 t^c = 1 < 1 , ;V/ ;J < 1 + r — = = 1 + 



Z X -^G) KE[F{a)\- E[F(a)] KE[F(&)} Z^G) 

The result follows from the extremely crude bound Z 1 - ^(G) > 8. □ 

We can now prove our main theorem. 

Theorem HI Let (3, 7 and A be efficiently- approximable reals satisfying ([I]). There is no 
FPRAS for DegreeFourPlanarTwoSpin(/3, 7, A) unless NP C BPP. 

Proof. We will give a randomised algorithm for PlanarCubicIS, using an FPRAS for 
DegreeFourPlanarTwoSpin(/3, 7, A) as an oracle (and also using the given FPRASes 
for /?, 7 and A). 

After receiving an instance G and h, our algorithm uses Lemma [25] which provides planar 
graphs J and J' with maximum degree at most 4 and also some approximation schemes for 
the reals K, 7 and A. With probability at least 1 — 1/15, these satisfy A > 4'^°-" and 
7 < \~\ V ( G *)\ < 1 and Equation ©. The algorithm then makes four calls to approximation 
schemes, suitably powered so that each call fails with probability at most 1/15. Thus, with 
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probability at least 2/3, Equation (0) is satisfied and all calls to the approximation schemes 
succeed. In that case, we will show how to determine (from the outputs of the approximation 
schemes) whether or not G has an independent set of size h. 
Let n = \V(G)\. By Equation ([6]), we have 

c -W K Z P,jA J ') < Z _ (G) < c l/4 g ^, 7 ,A(J') ^ 

Using the given approximation schemes for Zp a> x(J'), Zp rh \(J) and K, each with accuracy 
parameter 1/12 and failure probability 1/15, we can compute a value Z which , with probability 
at least 1 — 3/15, satisfies 

so 

e" 1 / 2 ^ < Z X - X (G) < e^Z. (14) 

Using the given approximation scheme for A with accuracy parameter l/2h and failure 
probability 1/15, we can compute a value A which, with probability at least 1 — 1/15, satisfies 
so 

e -l/2^ < yh < e l/2^_ (15) 

Suppose that all four calls to the approximation schemes succeed so that (fT4"J) and (|15[) 
hold. Recall that 

Z i,7,a( G )= E 7 c(<t) A £(ct) . 

u:V(G)^{0,1} 

If G has an independent set of size h then Z^- ^G) > X h so, plug ging in (fl~4"l) and ()15p . 

Z > e- 1 / 2 Z 1> _ >x (G) > e~ l / 2 \ h > e- 1 ^ 1 . 

Also, if G has no independent set of size h, then Z x - ^(G) < 2™ max(A h_1 , A n; y) < 2™A /l ~ 1 . 
So, plugging in (fill) and (fl~5j) and our lower bound for A, 



r\n\h 

Z < e^Z, - -JG) < e 1 / 2 2"A ?i - 1 < e V*±±- = e ^ 2 ~ n \ h < e2~ n \ h . 
— 1,7, A v / — — 4n — 

As long as n > 3, e _1 > e2 _n , so it is possible to determine from Z and A whether or not 
G has an independent set of size h. □ 

6 Approximating the log-partition function 

We start with a preliminary lemma, which will help us to show that our approximation is 
sufficiently accurate. 

Lemma 26. Suppose that f3, 7 and A are real numbers satisfying f3 > 1 > 7 > and A > 1. 
Then, for every planar graph G, Zp^ t \{G) > (1 + A)'^^'/ 4 . 

Proof. Let / be the largest colour class in a proper 4-colouring of G. Then / is an independent 
set of G of size at least |V(G)|/4. For every configuration a: V(G) — > {0, 1} which assigns 
spin to every vertex in V(G) \ I, w G {a) > A £ M. Thus Zp„^(G) > (1 + A)I 7 L □ 
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Our approximation algorithm is inspired by Baker's approximation schemes for optimisa- 
tion problems on planar graphs [I] . For a good explanation of her technique (which we use in 
our exposition here), see Borradailes's notes [6]. We will use the following notation (from [6]) 
to decompose a planar graph G = (V, E) which is embedded in the plane. We first define 
the level of each vertex. Vertices on the boundary of the embedding have level 0. Then, 
for i € {0, . . . , n}, the vertices with level i are those that are on the boundary on the graph 
formed from G by deleting all vertices whose level is less than i. For a fixed parameter k, and 
for every i £ {0, . . . , k — 1}, let V$ = {v £ V | The level of vertex v is equal to i modulo k}. 
Let Gi be the graph G — V{. By construction, Gi is {k — l)-outerplanar. Also, Bodlaender [2] 
had shown that every fc-outerplanar graph has treewidth at most 3k — 1. Also, this tree 
decomposition is easy to compute. Using a data structure of Lipton and Tarjan [20], Baker 
shows that the levels of vertices can be computed in 0(|U|) time. 

We can now prove Theorem [3l 

Theorem [3], Suppose that (3, 7 and A are efficiently approximable reals satisfying (3 > 1 > 
7 > and A > 1. There is a PRAS for PlanarLogTwoSpin(/3, 7, A). 

Proof. Consider input G = (V, E) with at least 3 vertices and an accuracy parameter e € 
(0, 1). Let n = \V\ and m = \E\ < 3n. Let f3 + , f3~ , X + , A - , 7 + and 7" be rational numbers 
(built into the algorithm) such that f3 + > (3 > f3~ > 1, 1 > 7 + > 7 > 7~ > 0, and 
A + > A > A~ > 1. Let k be any integer satisfying 

321og(2A + ) + 961og(/3+) 
£log(l + A-) 

Then let 

. 2nlog(2A+) 
= —— — . 

k[n + m) 

Using the given FPRASes for /3, 7 and A, compute /3, 7 and A satisfying e~ s (3 < (3 < e s (3, 
e _<5 A < A < e 5 A and e _5 7 < 7 < e <5 7. As in the proof of Lemma [25l adjust the output of the 
FPRASes to ensure 0+ > /3 > j3~ , 7+ > 7 > 7-, and A + > A > A". 

The first step is to compute a value Z satisfying 

Z < Z^-(G) < {2\ + f n/k {l3 + ) l2n/k Z. (16) 
This step is accomplished as follows. 

1. Using Baker's algorithm, construct the graphs Gi for i S {0, . . . , k — 1}. Each of these 
has treewidth at most 3(k — 1) — 1. 

2. Choose i 6 {0, ... ,k - 1} as follows. Let X = {i \ \Vi\ < 2n/k}. Note that \1\ > k/2. 
Now consider the 2m endpoints of edges in E. Choose i £ X so that Vi contains at most 
(2m)/\I\ of these. Note that \Vi\ < 2n/k and the number of edges with endpoints in 

is at most Am/k < \2njk. 

3. Use the algorithm of Yin and Zhang [271 Theorem 2.8] to compute Z = Zp~j(Gi). The 
running time of Yin and Zhang's algorithm is at most the product of a polynomial in n 
and an exponential function in the treewidth of Gi. In order to apply the algorithm, 
we first express the partition function Z^^j(Gi) as the solution to a Holant problem 
Holant(C/, J-) with regular symmetric T. See |27j for definitions and details. 
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4. Equation (|16|) now follows by noting that 

t: V(Gi)->{0,l} r':Vi-H0,l} 

where i(r') is the number of vertices u € V% with t'{u) = 1 and b(r, r') is the sum of 
the number of edges (n, with u S V^(Gj) and v € Vi and r(n) = t'(-u) = and the 
number of edges (u,v) with u £ Vi and u G V; and r'(u) = t'(d) = and c(r, r') is 
defined similarly (with spin 1). Then X^t'-v-^o 1} ^ T '/3 b ^ T ' T >y c ( T ' T > is at least 1 (since 

r can assign spin to every vertex in V^J and it is at most 2 ' (A) p 
To finish, note that 

e S(n+m) z ^ xiG) < Z MX (G) < e^+^^G), 

so since 

2n 

<5(n + m) < — log(2A + ) 
k 



and (from Lemma [ 

logZ^ )7)A (G) > (n/4)(l + A-), 

and 

(e/2)(n/4)log(l + A-) > y log(2A+) + — log(/3+), 



e- £ log(Z /3i7iA (G)) < log(Z /3>7jA (G))(l -e/2) 

2r? 1 2rf 

< log % 7>A (G) -5(n + m)- T log(2A+) - — log(/3 H 

< log Z^ A (G) - — log(2A+) - — log(/3 + ) 

< log Z. 



Similarly, 



logZ<logZ^(G) 

< 5{n + m) + log Zp nt \ (G) 
2n 

< T log(2A + ) + logZ ft7iA (G) 

< (l + e/2)logZ A7>A (G) 
<e £ logZ /3i7iA (G). 



□ 
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